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Abstract 

In this work we investigate characteristic modifications of the spectrum of cosmolog- 
ical perturbations and the spectral index due to chain infiation. We find two types of 
effects. First, modifications of the spectral index depending on interactions between 
radiation and the vacuum, and on features of the effective vacuum potential of the 
underlying fundamental theory. Second, a modulation of the spectrum signaling new 
physics due to bubble nucleation. This effect is similar to those of transplanckian 
physics. Measurements of such signatures could provide a wealth of information on 
the fundamental physics at the basis of inflation. 
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1 Introduction 



The current most popular models for inflation are based on chaotic inflation. In these 
models a scalar field rolls slowly subject to Hubble friction in a shallow potential. In 
[Ij we proposed an alternative scenaric0 that shares many of the features of slow roll 
chaotic inflation, but also differs in several important aspects. Our model is based 
on chain inflation and makes use of a series of several first order phase transitions. 
More precisely, we imagine a potential with a large number of small barriers that 
separate local, metastable minima. The barriers prevent the field from rolling, and 
without quantum mechanical tunneling the inflaton is stuck in a local minimum. 
By appropriately choosing the heights and widths of the barriers, one can obtain 
tunneling probabilities such that the field effectively rolls slowly down the potential 
through repeated tunneling events. In this way we can achieve a slow roll in the sense 
of having a slow change in ~ {p^ being the vacuum energy density), even if 
the potential for the fields is steep. The details of this process were worked out in 
[1], and it was also shown that suitable potentials might be find in flux compactified 
string theory. 

The main features of the model introduced in [1] is as follows. We assume that 
the bubbles, after their formation through tunneling, rapidly percolate and collide. 
The energy difference between two subsequent minima is temporarily stored in the 
bubble walls, and we assume that this energy is rapidly converted into radiation as 
the bubbles collide. In this way we obtain a coarse grained picture where the main 
effect of the barriers and the tunneling is to introduce a source term for radiation in 
the Friedmann equations. 



^For earlier attempts on models featuring inflation through chain of first order decays, see also 

[asiiiiii] 
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A scalar field can be understood as a fluid consisting of two components: a com- 
ponent corresponding to the kinetic energy, T, and a component corresponding to 
the potential energy, V. In the case of slow roll we have T ~ eV -C V, where e is 
the slow roll parameter, and as a consequence the dynamics is dominated by the po- 
tential energy leading to accelerated expansion and inflation. In our version of chain 
inflation the kinetic component is further suppressed relative to the potential energy. 
On the other hand, radiation is produced through the tunneling leading to prad ~ ^V- 
As a result we effectively have, to first order in e, a model consisting of a decaying 
cosmological constant and a coupled component of radiation. For the case of chaotic 
inflation, it is important to understand that it is the sub-dominant kinetic energy 
that determines the spectrum of the fluctuations. The kinetic energy corresponds to 
a hydrodynamical fluid with an effective speed of sound that is equal to the speed 
of light. In contrast, the potential energy does not correspond to a hydrodynamical 
fluid and lacks a well defined speed of sound. 

The amphtude of the primordial fluctuations are set by the speed of sound. The 
general result is 

where c^is the speed of sound of the hydrodynamical component. For chaotic inflation 
Cg = 1. In our model for chain inflation, the role of the kinetic energy is taken over 
by the radiation where = As a result, the primordial spectrum is corrected to 



p ^ (2) 

The result differs from chaotic inflation through a simple factor of v^- While this 
simple argument captures the main physics of the model, there are many important 
points of the derivation that are carefully discussed in p. 

In the present paper we discuss the possibility of further effects on the primordial 
spectrum from various sourcesH In the first part of the paper we will discuss the 
modifications to the spectrum of cosmological perturbations due to the presence of 
the non-negligible interaction between radiation and vacuum energy. We will discuss 
how they arise and appear to be specific to our model of chain inflation. In the second 
part of the paper we will instead consider how bubble nucleation affect the spectrum 
of perturbations, and in particular we will study the imprint of the size of the bubbles 
on the CMB. 



2 Effects due to interactions 

In [T] we derived a system of equations that determine the evolution of the comoving 
curvature perturbations during a period of chain inflation. The approach was based 

^Wc will leave aside the interesting possibility of having sizable contributions from isocurvature 
perturbations and/or non-gaussianities. 
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on the traditional analysis of scalar perturbations in field (slow-roll/chaotic) models, 
as presented in [B]. We start with a brief review of the approach used in pQ, and show 
that the method of [B] is not the most convenient one in the case of chain infiation. 
We will therefore propose another way of analyzing and re-writing the system of 
equations that is better suited for our model. 
We start from equation (97) in [1], 

• _ a(pjfp) 



(3) 



7^ = --J^—A- - aHirG^)^ - -HSy^ 



where 



a is the cosmological scale factor, H is the Hubble rate, is the gravitational 
potential, and G is Newton's constant 



• Qv is the energy- momentum transfer vectocl 

• a(f) = AixGHi, 

• TZ is the comoving curvature perturbation 

• p, p are the total energy and pressure density 

• e is the slow-roll parameter 

• is the comoving wavenumber for the perturbation 

• Sv,r is the relative entropy perturbation between vacuum (V^) and radiation (r). 

In the following we will neglect the term proportional to Sv,r- As a result our 
conclusions apply only to models with negligible contributions from isoentropic per- 
turbations, or, alternatively, just to the adiabatic component of the whole spectrum 
of perturbations. 

Comparing equations ([3]) with the analogous equations in [6] (in fiat space), we 
see the importance of interactions in our multicomponent system, as represented 
by the term —o^AttG^^. Note also that this can be conveniently re-written as 



■^See [T] for its exhaustive definition. For what we are concerned it is sufficient to define it through 
the equations 

p"" = Qv (4) 
= -AHp'- - Qv (5) 

where p^/^ is the energy density for vacuum/radiation. 
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—a'^AirG^ = a^H'^e. Following the literature (see for example [7|), it is customary 
to define a standard quantization variable <; = zTZ where 

The singularity at k'^ = a^H'^e is just an artifact of the choice of variables, as is 
evident from ([3])0. However, in order to better understand the implications for the 
spectrum of perturbations due to the new term, we choose to follow an alternative 
route using another change of variables. The (first order) action inferred from the 
equations of motion is given b}@ 

where O is a time-independent factor which, by comparison with known cases, was 
found to b^ O = k"^. By eliminating IZ in the action by using the first equation in 
(IH]), and defining 

= , 2 O (8) 

u = w^, (9) 

we obtain the action (up to a total derivative term) 

1 r 1 " 

S=- + -u(A + nh)u + —u^] dr] d\, (10) 

2 y 3 w 

and the equation of motion is 

///l,oloW;\ 

u+{-k^--'H^e u = 0. 11 

\3 3 w J 

Here we have used the conformal time defined as rj = J a~^dt and 7i = ^, where a 
prime represents a derivative with respect to the conformal time. An advantage with 
this way of writing the equations is the absence of any artificial singularities. 



2.1 The equation of motion 

// 

We will now expand — in slow-roll parameters such as e and the decay rates in such 
a way that we make sure to include contributions from derivatives of the slow-roll 

^This issue is absent in [6] in the case of flat space 

^This is the same action as in [B] with the derivative term integrated by parts. 
^In our previous paper we used O = —k^, but here we prefer use the other sign. 
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parameters up to order 0{e). We will study the general case of a vacuum energy e„ 
in the form of a power law, 

4 

en ~ (12) 
c! 

where the integer n labels the vacuum, and m/ is an energy scale (depending on cou- 
plings, extra-dimensions and similar features of the microscopical theoretical model). 
We will need the following formula (see appendix for derivation) 

£ ~ (^1 + He'^ - DF, (13) 

where DT depends on the decay rates per Tinit time T and is defined in the appendix. 
Note that we do not restrict to equal rates at every step. We find 



having neglected terms of order e'^^eDT, DT'^ and higher. Our equation of motion 
now becomes 



u" + Qa;^ - ^n'^e - n'^s - ^ (^(^1 +'^] - oHDV ) ) m = 0. (15) 
In a quasi-deSitter space, where 



■)) 



we have 



which allow us to read off 

The general solution for the equation 



reads 



u = ^ [c,{k)Hi'\-kri)+C2{k)Hi^\-kri)] , (20) 
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where hIP{x) and H^\x) are Hankel functions of the first and second kind, respec- 
tively. In the hmit — fcr/ ^ 1 we have that 



V —Tikr] V —T^kr] 

Following standard procedure we match the solution with the Bunch-Davies vacuu 



nndmg 



/2kc. 



u = ^.[^e'^'+^^^^y/^Hl^^-kT]). (22) 



For superhorizon scales {—kr] <ti 1) 



2_,5,,_ir(z/) 



so that, finally. 



Hl'^ ~ . -e-f 2^-i^(-A:r^)-^ (23) 



^^e-f(-i)2-fli^ (-^;^ ■ (24) 

r(f) 

This is what we need in order to obtain the spectrum of perturbations. 
2.2 Spectrum of perturbations 

The spectrum of perturbations is conveniently expressed through the comoving cur- 
vature perturbation, which is constant on superhorizon scales during inflation. We 
can obtain it using the first equation in ([3]), which we repeat here for convenience 

n = -T^A- (25) 



As a result we obtain 



H 1 / A; 

n= + (26) 



V2Mpiank£ y2F^ \aH 
and the spectrum becomes 



Plank 

From this we read off the spectral index 



2 f 2\ DV 

ns-l = l-2v ^ --e-2e-{l + -\e + — (28) 



^In the following section we will discuss this choice thorough, investigating the possibility of new 
fundamental physics at a scale larger than the Plank one. 
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which alternatively can be written as 



which is our final result. It is evident from this formula that the corrections to the 
spectral index due to the interactions between vacuum and radiation imply an extra 
tilt to the spectrum (blue or red depending on the value of DT). (Note that the result 
for D„V = Qy = is precisely the same as the usual chaotic/slow- roll result in the 
case of c = 2 as observed in [I].) It appears to us that this feature of the spectrum is 
strongly characteristic of a first order transition, since in deriving ([3]) in [Ij, we made 
use of specific aspects of first order transitions (such as the fact that the momentum 
perturbation for the vacuum was zero). 



3 Effects of a new intermediate scale 



3.1 The choice of vacuum 

In deriving our results for the spectrum of cosmological perturbations and the spectral 
index in the previous section, we have followed the standard procedure of matching 
our solution to the Bunch-Davies vacuum. Essentially this means that we resolve the 
issue of the non-uniqueness of the vacuum in a cosmological space time by tracking 
the modes to infinitely short scales, where the effect of cosmological scales such as the 
horizon can be ignored. At such scales there is a unique vacuum just as in Minkowsky 
space. This is the Bunch-Davies vacuum. 

As is well known there is a potential problem with this procedure since one can not 
reliably track the modes to scales shorter than the Planck (or string) scale without 
taking into account effects of string theory and quantum gravit}|§. Hence there are 
likely corrections to the choice of vacuum of order ^, where A is the scale of new 
physics. This is known as the transplanckian problem. Actually, it represents more 
of an opportunity than a problem since it could be an observational window to new 
physics. 

In our case we have yet another scale that enters. We have assumed a coarse 
graining over the nucleating bubbles that is valid only for scales substantially larger 
than the size of the bubbles, r^. Hence, we have full control over the evolution of the 
perturbations only while their wavelength is larger than the size of the bubbles. If 
we follow the evolution of a specific mode backwards in time it will eventually reach 
a scale as short as the size of the bubbles, and our picture breaks down. 

What is the effective quantum state that should be used as an initial condition 
at this point? It is in general very difficult to give a precise answer to this question 
both because of the usual difficulties due to quantization in curved spaces, and also 

«See for example [H [9l [IHl El [l2] 
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because of the great generality of our model of chain inflation (no field theory model 
is specified). Without a detailed model, we have no other option than to impose an 
effective initial condition for the perturbations, and to postulate their creation out of 
the vacuum. This is formally very similar to the case of the transplanckian problem. 
In several works (O [9], [101 [HI 112]) it has been argued that initial conditions must 
be imposed at the Planck scale (or string scale) due to our ignorance of physics at 
higher energies. In our model the scale for new physics will instead be the size of the 
bubbles, but the analysis, that we now review, will be more or less the same. 

We begin by noting that the physical momentum p and the comoving momentum 
k are related through 

^ = «P=-;^' (30) 

where rj is conformal time, p is the physical momentum, k is the comoving momentum 
and a is the scale factor. We impose the initial conditions when p = A, where A is the 
energy scale important for the new physics given by A ~ l/r;,. In our case Vb is just 
the size of the bubbles. We find that the conformal time when the initial condition 
is imposed to be 

"» = -m (31' 

As we see, different modes will be created at different times, with a smaller linear size 
of the mode (larger k) implying a later time. 

In our case we would in principle be able to calculate the form of the perturbations 
at the the scale r^, by tracing the evolution backwards in time, through the nucleating 
bubbles, to even smaller scales. Presumably the result would depend on the fine 
details of the physics of bubble nucleation, which is beyond the scope of the present 
paper. 

Instead we will take the same attitude as in [10] and encode the unknown new 
physics into the choice of the vacuum. The claim of (TUj is that the primordial 
spectrum is corrected through a modulating factor. These results can be directly 
taken over to our case with the result that 

where = ^ is the speed of sound. In the transplanckian case, A is typically constant 
and equal to the Planck scale or string scale. The modulation of the spectrum comes 
from the rolling inflaton that leads to a changing H. In our case, A will also be 
changing, but the amplitude ^ is nevertheless expected to be small, and not to 
change very much during inflation. The argument of the sine, on the other hand, is 
a large number and can easily change by several times 27r during the relevant time 
period for the generation of the primordial perturbations. 

Let us now investigate in more detail what the effect will be in the case of chain 
inflation. 
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3.2 The size of the bubbles 



To proceed we need to know more about the process of nucleation of bubbles. The 
rate of bubble formation per unit time and physical volume, in the analysis of |14j . 
[15], is given by the exponential of the euclidean instanton action, Se, responsible for 
the tunneling (the so-called "bounce"), 

r-e"*^. (33) 

The action evaluated on the bounce is given bj{^: 

5*^ = -yr^Ae + 27rV'^, (34) 

where r is the radius of the bubble, Ae the change in energy due to the nucleation 
of the bubble, and 5* is the bubble's wall tensionj^i In principle, there is also a third 
term present due to the effect of gravity, but we assume the size of the bubbles to 
be much smaller than the Hubble scale so that we can ignore it. The critical radius 
that allows for the nucleation of a bubble that will successfully expand, and therefore 
enables tunneling, is obtained by extremizing the above Euclidean action. The result 
is 

n = '"critical = ('^5) 

Se\ = 5-- (36) 

'^lr-=rcritical 2 (Ae) 

The setup outlined in these formulas is a static one. The tunneling occurs between 
two vacua of the theory, and any possible time evolution of the background is not 
taken into account. In our scenario, on the other hand, we have a chain of tunneling 
events occurring through time. The length scale signaling new physics (corresponding 
to the radius of the nucleated bubbles) depends on the time when the particular mode 
of interest is produced. For simplicity, however, we will assume that the change in 
the radius is slow enough that we can use the above analysis. 

Accounting for the time evolution of the background, when computing the critical 
radius of the bubbles at a given time, is most easily achieved by expressing the 
variation of the energy density due to the nucleation as 

Ae~-^(r) = 6i/3M|,,,,(r)5, (37) 

where we have used the Friedman equation, and defined (r) = (FV^^ to be the average 
tunneling time (we recall that F is the decay rate per unit time)LJ. Also, the surface 

^We recognize in this the variation of the Gibbs energy for the nucleation of a bubble: with a 
different normalization for what concerns energy density and surface tension. 

^"if the tension is due to a scalar field, (j), we have that S — J d(j)VV, where V is the potential. 
^^All averages are taken with the distribution of vacuum phases p^^ — emPm{t), see appendix and 

m 
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tension 5* needs to have a time dependence. It is more convenient, though, to express 
S through the extremized action, as 



Eventually we find 
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cvr^ / n \Mpi 



where 7] ~ 2.5 ■ 10~^ from the normalization of the spectrum. With e = 10~^ we find 

nH ~ 3.9 • IQ-^ Slc-Ur)-U-)-^. (40) 

n 

In our calculation we have ignored time dependent corrections to, e.g., Se that are 
suppressed by 1/n. 

Let us now consider the possible observational implications of the above effects. 
Successful chain inflation requires that while the vacuum undergoes the transitions, 
the phase distribution in the universe is peaked consecutively on the various phases. 
That is, the transitions occur consecutively, and in a short time (shorter than the 
Hubble time). Rapid tunneling implies that Se should be at most of order one, in 
order for {t)H <^ 1. If we then use e = ^{^) = 2{^{'^)(^) ^ we find that (t)(^) 
needs to be small. With at peaked distribution we have, to a good approximation, 
(t)(^) ~ and we see that n needs to be large. 

Turning back to equation fl32|) . and the corrections to the spectrum from the 
presence of a new scale, we know from the work on transplanckian physics that values 
of the order e ~ 10~^ and ^ ~ 10~^ could possibly yield an observational effect. The 
restriction on e comes from the requirement that H changes in an appreciable way in 
order for there to be a modulation. Using H ~ and fl5^ we have, following 



Ak ttH 

(41) 



k eA 
hat ^ of 

k 

on the order of percent are easily obtainable within our model using values of n ~ 10 



where, in our case, A ~ l/r;,. We see that ^ of order one, and a reasonable amplitude 



4 Discussion 



As we have argued, chain inflation will lead to several new effects on the spectrum of 
primordial perturbations. In particular, our calculations show that the spectral index 
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is changed from the naive one due to the presence of interactions between radiation 
and the vacuum energy (the contribution proportional to Qy oc e in formulas fl28ll29p ). 
The detailed predictions depend in a sensitive way on the distribution of the decay 
rates between the different minima. If it would be possible to measure these decay 
rates in a precise way, they would provide us with a wealth of information on features 
of the (effective) potential of the underlying fundamental theory. One needs to keep 
in mind, though, that it is necessary to distinguish these effects from other similar 
effects that could arise from non-standard potentials in other models of inflation. 

Another, possibly more characteristic prediction, is the existence of signatures 
similar to those that could be generated through transplanckian physics. That is, 
a modulation of the spectrum due to the presence of a fundamental scale. In case 
of transplanckian physics, it is the Planck scale (or string scale) that determines 
the effect, while in the case of chain inflation it is instead the size of the nucleating 
bubbles. It is interesting to note that the model quite naturally, without much fine 
tuning, gives rise to effects of a reasonable magnitude that possibly could be detected. 
As in the case of transplanckian physics we have only been able to make a very rough 
estimate of the size of the effect. 

In order to make better predictions of observational signatures, a precise model 
with an explicit potential and field content needs to be specified. In [Ij, based on work 
in [16] and [17], we proposed that flux compactified type IIB string theory provides 
such models in a natural way. In that work we focused on the stabilization of the 
complex structure moduli using fluxesH With the help of monodromy transforma- 
tions, generated by going around singular points in the moduli space of Calabi-Yau 
compactifications, we were able to show the existence of long sequences of minima 
of the necessary form. A quadratic behaviour, with c = 2, typically arises when the 
axiodilaton is stabilized independently of the complex structure moduli, while the 
linear behaviour with c = 1 arises when the axiodilation is stabilized together with 
the complex structure moduli. The detailed form of the potentials depends heavily 
on the choice of Calabi-Yau manifolds and fluxes, but the overall features seem to be 
rather generic. In particular, the barriers in between the minima are expected to be 
such that an effective slow roll behaviour arises. 

It would be interesting to further explore the possibility of generating potentials 
for chain inflation through string theory. Given the difficulty in finding appropriate 
potentials for standard inflation, we believe this to be a worthwhile enterprise. 



In our simplified model the Kahlcr moduli, i.e. the moduli determining the sizes of the extra 
dimensions, were assumed to be fixed by other physics. 
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A The slow-roll parameter and its first derivative 
in a power-law chain inflation model 



During inflation we naturally expect: 



anc^^' 



~ (42) 



^ y^emPm- (43) 



Here Pmif) is the fraction of volume occupied by the vacuum m at time t and its time 
evolution is given by (see [T]) 

prn = -TmPm + Tm+lPm+1- (44) 



From this, we find 



and from (fT2l) we see 



^ AemTmPm, (45) 



e 



Ae„ ~ c— (46) 



m. 



Then, using this and f j42l) . we find for the slow- roll parameter 



If we now define an average (^) as 



En ^nPn 



the slow-roll parameter is given by 



2ii n 



The time derivative of the slow-roll parameter is as follows. From (H71) and using ( l46l) 

'"2^n^+2ifl En^nPn (En^nPny ' ^^^^ 



^•^Here and in the following represents the energy density of the vacuum in the interior of the 
bubbles. We in general expect this formula for the following reason. First of all during inflation 
the total energy density is dominated by the vacuum component. The latter is then given by the 
contributions respectively of the interior of the bubbles and the walls. But the energy density of 
uncoUided walls is proportional to the energy difference between two consecutive vacua, while the 
one of the interior of bubbles is proportional to the energy level, which is greater than the difference. 
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Let us focus on the two terms in each of the brackets. The second one is simply 
It is the easy to see that using the numerator in the first term reads 

^ ^m^mPm ^ ^ = ^ m+lPm+1 ( Tm —r'^m+l ) 

^-^ ^ — ^ Km m + 1 / 

= ~ rm+lPm+1 ( ^ ( 7— r ) H A(rm+: 

V \ m + 1 / m 



We find that: 



If we now define 



^ , 2\ / ~ \ 2 

o / , r\ \ /r 



Then, from ( jSOil and (H9l) . we have 



(52) 



where we have defined, for any quantity 

^ (/m) = fm fm—l ■ (^3) 



m + 1/ (m + 1)^ 



/ 9\ r r ~ r r af f 

8^(1 + -) He'-{c-l)eaU-)-'-e{-AT){-)-' + {c-l)e{ -){-y\ (56) 

\ cj n n n n nn — 1 n 



where all the averaging has been made using the distribution = e^P 
For ease of notation, we define 



m- 



~ F F ~ F F AF F 

DF = (c- l)ea2f (-)-i + .(-AF)(-)-i - (c- 1)5( -)(-)"\ (57) 

^rn n n nn — I n 
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